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Abstract We present the idea of a Mott- Anderson freeze-out that suggests a key role of the localization 
of the hadron wave functions when traversing the hadronization transition. The extension of hadron 
•/^ ' wave functions in dense matter is governed by the behavior of the chiral quark condensate such that 

its melting at finite temperatures and chemical potentials entails an increase of the size of hadrons 
and thus their geometrical strong interaction cross sections. It is demonstrated within a schematic 
(~| , resonance gas model, that a kinetic freeze-out condition reveals a correlation with the reduction of the 

^ i' chiral condensate in the phase diagram up to 50% of its vacuum value. Generalizing the description of 

the chiral condensate by taking into account a full hadron resonance gas such correlation gets distorted. 
We discuss, that this may be due to our approximations in calculating the chiral condensate which 
I disregard both, in-medium effects on hadron masses and hadron-hadron interactions. The latter, in 

particular due to quark exchange reactions, could lead to a delocalization of the hadron wave functions 
in accordance with the picture of a Mott- Anderson transition. 

Keywords Chiral condensate • chemical freeze-out • Mott- Anderson transition 
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lO ' 1 Introduction 

' The investigation of the structure of the QCD phase-diagram is one of the goals of heavy-ion collision 

experiments as well as theoretical studies within lattice QCD (LQCD) and effective field theory ap- 
proaches to the non-perturbative sector of QCD. Of particular interest are the conditions under which 
the approximate chiral symmetry of the QCD Lagrangian get restored and whether this transition is 
. to be necessarily accompanied by the deconfinement of quarks and gluons. More detailed questions to 

' the QCD phase transitions concern their order, their critical exponents as well as the existence of a 
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critical point or even a triple point in the phase diagram. Most promising tools for the experimental 
determination of these characteristics are the energy scan programs at CERN, RHIC and at the up- 
coming dedicated facilities of the third generation: FAIR and NICA. The systematic analysis of higher 
moments of distributions of produced particles in their dependence on the collision energy and the size 
of colliding systems shall provide answers to the above questions and allow for direct comparison with 
predictions from the underlying theory, as provided by LQCD, see, e.g., Ref. 

As long as the applicability of LQCD methods is bound to the region of finite T and fis/T ^ 1, any 
predictions for the phase structure of QCD at high baryon densities including the possible existence 
of critical points will rely on effective models. To be relevant for the discussion of the above problems 
these models have to share with QCD the property of chiral symmetry and its dynamical breaking as 
well as a mechanism for confinement and deconfinement. 

At present, the position of (pseudo-)critical lines in the QCD phase diagram is constrained by the 
chemical freeze-out parameters {T^ , jig) which have been obtained from the statistical model analysis 
of particle yields obtained in heavy ion collisions 0; Q . One of the most striking observations is the 
systematic behaviour of these parameters with collision energy 0; [H; 1 which has recently been 
given a simple parametric form [7[ . It has been observed, that the resulting freeze-out curve in the phase 
diagram is closely correlated to the thermodynamical quantities of the hadron resonance gas described 
by the statistical model, see Fig. [TJ These phenomenological freeze-out conditions make statements 
about the mean energy per hadron {E)/{N) ~ 1.0 GeV, the dimensionlcss entropy density s/T^ ~ 7 
and a total baryon and antibaryon density hb + rig ~ 0.12 fm"''. The freeze-out line provides a lower 
bound for the chiral restoration and deconfinement transitions in the phase diagram. Being coincident 
at low densities as inferred from LQCD Q, both transitions need not to occur simultaneously at high 
densities, thus allowing for an island of a quarkyonic phase between hadron gas and quark-gluon 
plasma with a (pseudo-)triple point [loj . The question appears for the physical mechanism which 
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Fig. 1 Freeze-out parameters from the statistical model [T3| compared to phenomenological rules for explaining 
their systematics in the temperature - chemical potential plane. Ellipses indicate regions where dedicated future 
experiments (FAIR-CBM and JINR-NICA) could provide new data for more detailed investigations. 



governs the chemical freeze-out and which determines quantitatively the freeze-out parameters. One 
aspect is provided by the requirement that hadrons should overlap in order to facilitate flavor exchange 
reactions which establish chemical equilibrium. This geometrical picture of freeze-out is successfully 
realized in a percolation theory approach [llj. Another aspect is the dynamical one: when the equation 
of state (EoS) possess softest points (e.g., due to the dissociation of hadrons into their quark and gluon 
constituents with a suflacient release of binding energy involved) which is correlated with the freeze-out 
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curve in the phase diagram, then it is expected that the hadron abundances are characterized by the 
corresponding T and /is values [l^. The dynamical system got quasi trapped at the softest points 
for sufficient time to achieve chemical equilibration before evaporating as a gas of hadron resonances 
freely streaming to the particle detectors. Also the kinetic aspect of fast chemical equilibration was 
discussed in the hadronic gas when accounting for a multi- hadron dynamics . 

All these mechanisms are appealing since they provide an intuitively clear picture, but they are 
flawed by the fact that their relation to fundamental aspects of the QCD phase transition like the 
chiral condensate as an order parameter are not an element of the description. 

In the present contribution we develop an approach which relates the geometrical as well as the 
hydrodynamic and kinetic aspects of chemical freeze-out to the medium dep endence of the chiral 
condensate. Wc demonstrate within a beyond- meanfield extension p^ : [Tsl : [l6l : It of the Polyakov 

NJL model [20,; 21; 22; 23J how the excitation of hadronic resonances initiates the melting of the chiral 
condensate which entails a Mott- Anderson type delocalization [l^; [H; [1^ of hadron wave functions at 
the hadronic to quark matter phase transition [23]. Already for a schematic resonance gas consisting 
of pions and nucleons with artificially enhanced numbers of degrees of freedom we can demonstrate 
that a kinetic freeze-out condition for the above model provides quantitative agreement with the 
phenomenological freeze-out curve. 

We propose a generalization of this approach by including a more complete description of the chiral 
condensate which accounts for contributions from the hadron resonance gas and the chiral dynamics 
implemented from the PNJL model. We demonstrate, that in this approach correlations between freeze- 
out and chiral condensate get distorted. We discuss, that this may be due to our approximations which 
exclude the in-medium effects on hadron masses as well as the hadron-hadron interactions which could 
result in delocalization of the hadron wave functions in accordance with the picture of a Mott- Anderson 
transition. 



2 Mott- Anderson freeze-out and chiral condensate 

Recently, we have sug gested, that the chemical freeze-out in heavy ion collisions can be linked to the 
chiral condensate |28|. The starting point is based on the assumption, that the chemical freeze-out 
sets in if the collision and the expansion time arc equal. Thus, in the temperature-chemical potential 
(T, /i)-plane, the kinetic freeze-out condition is formulated as 

Toxp(I', At) = Tcoii(r, ^J) , (1) 

where Tcxp(T, ^) is the expansion time of the hadronic fireball, and 

'^cMiT,^^) = X^cr^jJ^j , (2) 

i-.j 

is the inverse relaxation time for reactive collisions with i,j — tt, N, ... running over all particle species 
in the hadron resonance gas. For the hadron-hadron cross sections aij we adopt the geometrical Povh- 
Hlifncr law, [H; [s^ 

cT.,^X<jj){r^) , (3) 

where A is a constant, being of the order of the string tension, A ~ 1 GeV/fm = 5 fm^^. Note, that 
this behavior was also obtained for the quark exchange contribution to hadron-hadron cross sections 

. . .. 

A key point of our approach is that the radii of hadrons depend on T and fi and that they diverge, 

when hadron dissociation (Mott-Anderson delocalization) driven basically by the restoration of the 
chiral symmetry, sets in. This has been quantitatively studied for the pion [3^ . where it has been 
shown, that close to the Mott-Anderson transition the pion radius is well approximated by 

rliT,f,) = ^F-'{T,^^) . (4) 

In addition, it has been demonstrated, that the GMOR relation holds out to the chiral phase tran- 
sition, where pions merge the continuum of unbound quark matter (33l : [3^ . Furthermore, since the 
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current-quark mass is T— and /i— independent and the pion mass is "chirally protected" , the T— and 
/^—dependence of the chiral condensate should be similar to that of the pion decay constant, 

Fl{T,^i)^-m^{qq)T^Jml. (5) 

Consequently, the relationship between the pion radius and the chiral condensate reads 

'^^(T,M) = ^I(^<7>T,,|-^ (6) 

The Mott- Anderson delocalization of the pion wave function due to melting of the chiral condensate, 
as expressed in this formula, is the important element of the hadronic freeze-out mechanism suggested 
in [13. 

For the nucleon, we assume that its radius consist of two components; a medium independent hard 
core radius ro, and a pion cloud contribution as 

rl{T,^i)=rl+rl{T,^x) , (7) 

where from the vacuum values = 0.59 fm and = 0.74 fm one gets tq = 0.45 fm. 

For the expansion time scale we adopt the relation which follows from the entropy conservation, 
S = s(T, /i) y(Toxp) = const with V{tcxp) being a fireball volume. Assuming, that y(Tcxp) oc r^^p one 
gets 

Texp(r,M) = a s-i/3(T,/i) , (8) 

with a being a constant of the order one. 

To quantify the freeze-out conditions in the {T — fi) plane, we first consider a medium which is 
composed of pions and nuclcons only, however, with a variable number of degrees of freedom. In this 
way, we account effectively for additional mesonic and baryonic stats in the hadron resonance gas. A 
generalization to a medium composed of all hadrons with physical degrees of freedom will be discussed 
later. 

We start from a medium-modified chiral condensate obtained within the chiral perturbation theory 

(gg)vac mlF^ ' 

where the pion- nucleon sigma-term is S-^n ~ mo(9mjv/3mo) ~ 45 MeV [S^], and the scalar densities 
of pions Us^TT and nucleons Us^n read 

ns,N^^l^^dpp'^^^{f+iEN{p)) + f],{ENip))} ■ (11) 

with (Itt the pionic and the nucleonic number of degrees of freedom. In the chiral limit and for (1,^ — 3 
we have, ns.Tr/mTr = T'^/8. Fig. [5] shows the chemical freeze-out lines in the (T — /x)-plane obtained 
from the kinetic freeze-out conditions ([1]) for three sets of pion and nucleon degrees of freedom. The 
best description of the phenomenological freeze-out points is obtained for = 8 and djv = 20. The 
inset of Fig. [5] shows a dependence of the position of the freeze-out line on the expansion (solid-line) 
and the collision (dashed-line) times used in the model for ~ 8 and d^ ~ 20. 

From Fig. [5] one concludes, that the chemical freeze-out in the energy range between RHIC and 
GSI (SIS) occurs 3-5 fm/c after the expansion of the collision fireball has started. It is also evident 
from Fig. [2l that increasing d^^ results in decreasing freeze-out temperatures in the meson-dominated 
region at small fis/T. On the other hand, increasing dj^ shifts the entire freeze-out curve towards lower 
T and fiB- The freeze-out curves calculated in this simplified model correspond to a reduction of the 
chiral condensate up to 50% of its vacuum value, c.f. the upper right panel of Fig. [3] Thus, according 
to Eq. dS]), the hadron-hadron cross section at freeze-out in ([2]) is roughly 3-4 times its vacuum value. 

In the next section, we introduce the first results of extending the above freeze-out model including 
correctly all hadrons contained in the hadron resonance gas. We start our discussion by introducing a 
generalization of the chiral condensate. 
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Fig. 2 Freeze-out curves according to the kinetic freeze-out condition, Eq. for the schematic resonance gas 
model with different number of pion (d^) and nucleon (djv) degrees of freedom compared to phenomenological 
values (symbols) The lines correspond to different choices for d,r and cIn given in the legend. The inset 
shows the lines of constant expansion time (black solid) and constant collision time (red dashed) for values 2, 
3, 4, and 5 fm/c. Lines for Tcoh are so close together that their behavior determines the location of hadronic 
freeze- out. 




3 Chiral condensate in a hadron resonance gas 

The chiral condensate is introduced based on the thermodynamical potential f2{T,fi) which is decom- 
posed into contributions from quarks and gluon (mean-field) and the hadronic (quantized) fiuctuations 
in the channels of Goldstone bosons (G = tt, K, t], 77'), mesons {M = /o, p, lo, ...) and baryons {B = 
N, A, S, A •••), 

^{T,^l) = QMF{T,fi) + r2G(r,^) + Qm{t,h) + QB{T,fi) . (12) 

For the mean-field contribution in the quark and gluon sector we employ the expression obtained 
within the PNJL modeQ [l^ , 

I2mf(T,/.) = -2^11^ {3s f + T In [N+isf)N^ief)] } + ^ - ^| +um, (13) 



where N^{e) ^ 1 + S^e"'^'-^^''^ + 3<?e-2/9(eTA') + g-3/j(eTM)^ ^ith e/(p) = ^ + mj {T , fi) being the 

quasi-particle energy and •mf{T,ii) = m^j + af{T,n) the dynamically generated mass of quark with 
flavor /. The vector fields ujf renormalize the quark chemical potentials Hf = jif^ — Uf. For the 
Polyakov-loop potential we take the logarithmic form motivated by the SU(3) Haar measure, 



- i a{T) -P^ + b{T) ln(l -d-P^ +8<P^ -3 <P^) 



, (14) 



with the corresponding parameters a{T) and b{T) from Ref. [ig. A possible dependence of the Tq on 
the number of active flavors and on the chemical potential has been suggested in Ref. [3]. Here we 
use the constant value of To = 187 MeV appropriate for the (2-|-l)-flavor model. 

In the following we consider only the on-mass-shell hadron contributions and neglect the continuum 
correlations. This approximation is motivated by the fact that hadronic freeze-out occurs in a region 
of the phase diagram which docs not exceed the limits of the hadronic phase. It may, however, turn 



^ We use the approximation $ = $ = [1 + 2 cos(03/T)]/3 for the mean value of the traced Polyakov loop. 
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TT 


92.4 


3 


1 


K 


113 


4 


1/2 


V 


124 


1 


1/3 


v' 


107 


1 


2/3 



Table 1 Parameter values for Goldstone bosons, see also [44j |. 



out to be insufficient in tlie vicinity of tlie chiral restoration transition where the energy gap between 
hadronic bound states and the threshold for the quark continuum states becomes too small to suppress 
the excitation of the latter relative to the former. 

The Goldstone and meson contributions to the thermodynamical potential are given as 



G,M 



E 

G = ir,.. 

M=fo,. 



(27r)3 



G,]\I 



(fc) 



T In 



1 



(15) 



where the index G (A/) denotes the actual meson with degeneracy factor da {du) and the dispersion 



fc2 



The contribution of baryons to the partition function is considered as a quasiparticle Fermi gas 

%{EB{k)+Tl. 



B=N,... 



(2t 



}}-(16) 



In this work we study flavor-symmetric matter where quantities for the sector of the light flavors 
u and d are denoted by the subscript q. We define, aq = (tT„ + ad)/^, while m = rriu = md and 
rriQ = rriQ^u — THo^d- The light quark condensate is defined as 



dd) = 



d 



f2{T, 



(17) 



In this work we will not yet solve selfconsistently the gap equations for the set of order parameters 
{tT/}, {ujf}, ... which correspond to the thermodynamical equilibrium state characterized by the global 
minimum of the thermodynamical potential (|12p . i.e. including contributions from hadronic resonances. 
Such contributions have been discussed for mesonic fluctuations in [s^; [40| for the NJL model and in 
[ssl : [4ll | for the PNJL model. Also recently mesonic fluctuations were included in these models within 
the functional renormalization group approach [42|; 111] . 

To get dniM.B / drriQ j we use rather simplified mass formulas 

tob = (3-iVs)(o-g + mo) + A^s(o-s+?7io,s) + KB , mA/ = (2-Afs)((Tq + mo)-l-Afs(o-6+mo,s) + KA/ , (18) 

where kb and km are state-dependent constants for baryons and mesons, respectively. 

Let us first discuss a derivation of the two-quark condensate from Eq. (jl7p including contributions, 
beyond the quark mean-field, of mesons and baryons forming the hadron resonance gas. To this end, we 
employ the schematic mass formulas for baryons and mesons p^ . the GMOR motivated relationship 



d 



d 



dniQ u dm, 



0,d 



and assume flavor symmetry of the vacuum (ss) 
and Ref . H . We obtain 



-47V, 



dpp^ m 



[i-fi-u] 



dd)/2. For the coefficients tg, see Table [1] 



2^2 2cgF^Eg{p) 



^ + fG{EG{p)) 



E 

M=fo,ui,.. 

E 

B=N,A,... 



dM(2-iVf) 
dB(3-iVf) 



dpp'^ ruM 
2^2 Em{p) 



^ + fM{EM{p)) 



dpp^ 



rriB 



27r2 EBip) 



[l-f+{EBip))-fB{EBip))] 



(20) 
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where cq = Fq / and the PN JL model distribution functions for quarks read 

/±(£) = {,^[e-/3(^'T/^) + 2e-2^('^=FM)] + g-3/?(eTM)}/7v±(e) . (21) 

The values of the Polyakov loop ^ in the phase diagram are found from a similar gap equation cor- 
responding to the solution of the extremum condition dH/d'P — 0. The distribution functions for 
mesons and Goldstone bosons as bound states of quark with flavor /i and antiquark with flavor /2 
are: Jm.g — ^/ {'^^v{P{E'M,g ~ ~ M/2))] ~ !}• The distribution function for baryons composed of 
three quarks with flavors /i, /2 and /s are: = l/{exp[/3(i?B T /^/i T M/2 T M/3)] + !}■ 

We will use now the property of chiral protection of the pion mass in the phase with broken chiral 
symmetry which entails that the in- medium modification of the chiral condensate shall be compensated 
by that of the squared pion decay constant 

(gg) _ (gg)vac _ _rr^ 

~ F^ ~ mn ' ^ ' 

In the following we drop the subscript index "vac" and use F^^ for the vacuum value of the pion decay 
constant. The final form for the in-medium modification of the chiral condensate is then 



■ E 



'^G^'G f ^p^^f^(^E^(^p^y (23) 



An^F^J Eg{p) 



The above result for the two-quark condensate reproduces that found in Ref. [44| for the hadronic 
resonance contributions, and generalizes it by including the quark mean-field contribution according 
to the chiral quark dynamics from the PNJL model. 



4 Results of the PNJL-hadron resonance gas model 

In the following, we summarize first results for the chiral condensate, which is the main input to 
our description of the freeze-out conditions within the Mott- Anderson hadron localization mechanism. 
For these exploratory studies, we neglect the hadron effective mass and the rcnormalization of the 
chemical potentials, so that we stay rather close to the standard statistical model approach which 
operates without medium modification of the hadronic states. 

We use standard values of the NJL model parametrization, e.g., from (4^, with GA^ = 2.316, 
A = 602.3 MeV, toq = 5.5 MeV, ms,o = 138.7 MeV, = 140 MeV, Mk = 495 MeV, = 92.4 MeV 
and fK ^ 93.6 MeV. 

In the first step, we confirm consistency with our model assumption, that quark degrees of freedom 
do not yet interfere at the freeze-out temperatures and chemical potentials. In order to do that we 
neglect in the in-medium condensate (j23p contributions from Goldstones as well as mesons and baryons 
by setting the corresponding degrees of freedom to zero, du = ds = da = 0- The remaining expression 
is that for the PNJL model condensate in the mean-field approximation, 

7^-1-5^ - -.,-4ivJf^^ [/+ + /,] . (24) 

The results for the freeze-out lines with the above condensate are shown in the upper-left plot of 
Fig. [21 They confirm, that the reduction of the vacuum condensate is still by less than 10% along 
the chemical freeze-out data. Since the modification of the mean-field contribution stemming from the 
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Fig. 3 Chiral condensate in the T-jj, plane compared to freeze-out data. Upper left: two-quark condensate for 
PNJL quark mean fiel d without hadronic contributions; upper right: schematic resonance gas model [231; lower 
left: present work (|23p excluding quark meanfield [4j including hadronic states with masses up to 1 GeV; lower 
right: the same for masses up to 2 GeV. 



quark excitations is proportional to /<j, it is effectively suppressed by the Polyakov-loop, compared to 
the standard NJL model case. 

Wc investigate further the role of the hadronic states for the behavior of the chiral condensate. First 
we consider a schematic model by reducing a complete set of states to nucleons and pions only. Their 
numbers of degrees of freedom, = 8 and cIn — 14 are chosen such that to get the best description of 
the freeze-out data. We observe that the freeze-out curve is correlated with the condensate reduction 
to 50 %. 

The lower panels of Fig.[3]show results for the quark condensate from Eq. when the contribution 
of the quark mean-field is excluded. This assumption corresponds to the two-quark condensate form 



4J], including hadronic states with masses up to 1 GeV (lower left) and up to 2 GeV (lower right). 



We observe that the shape of the lines of constant quark condensate in the (T- /i)- plane is different 
from the systematics of the freeze-out curve and thus different from the behavior of the schematic 
resonance gas model. Extending the mass limit of hadronic states from 1 GeV to 2 GeV has only a 
minor influence on the onset of the chiral restoration (reduction up to 50%), visible in the region of 
large baryon-chemical potentials, since in this mass range predominantly baryonic states are added. 
Obviously, the straightforward generalization of the simple kinetic freeze-out model based on a scaling 
of hadron cross sections with the chiral condensate towards a realistic resonance gas model requires 
improvements of the assumptions made in the present work. We will discuss them in the concluding 
section. 
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5 Conclusion and outlook 

We have developed a chemical freeze-out mechanism which accounts for a strong medium dependence 
of the rates for inelastic flavor-equilibrating collisions. Wc based our concept on the delocalization of 
hadronic wave functions and growing hadronic radii when approaching the chiral restoration. This 
approach relates the geometrical (percolation) as well as the hydrodynamic (softest point) and kinetic 
(quark exchange) aspects of chemical freeze-out to the medium dependence of the chiral condensate. For 
our model calculations we have employed a beyond-mean-field, effective extension of the Polyakov-loop 
NJL model. We could demonstrate how the excitation of hadronic resonances initiates the melting of 
the chiral condensate which entails a Mott- Anderson type delocalization of the hadron wave-functions 
and a sudden drop in the relaxation time for flavor equilibration. 

For a schematic resonance gas consisting of pions and nucleons with variable numbers of degrees of 
freedom we could demonstrate, that the kinetic freeze-out condition can provide quantitative agreement 
with the phenomenological freeze-out curve for suitably chosen values of the degeneracy of pionic and 
nucleonic states. 

We have developed this approach further by replacing the schematic resonance gas model for the 
chiral condensate with a full treatment of hadronic resonances and studied their influence on the chiral 
condensate. To this end generalized models for hadron resonance masses were required, which reveal 



their dependence on the current quark mass like, e.g., the model of Leupold [44| which we have adopted 
here. 

The obtained result for the chiral condensate which triggers the freeze-out in the Mott- Anderson 
localization model, showed that the approach by Leupold to the chiral condensate in a hadron resonance 
gas, does not show the expected correlation between chiral condensate and freeze-out. 

We see the main reason for the failure of the improved model in the inconsistency of the Mott- 
Anderson delocalization picture with the assumption of on-shell hadronic spectral functions in the 
thermodynamical potential (|T^ as well as in the very crude baryon and meson mass formulae ((T5|) . 
The latter provides a value for the pion-nucleon-sigma term which is about a factor three too low. 
Clearly, the effects of chiral hadron-hadron interactions arc missing in (IT^ . At high phase-space oc- 
cupation, one expects, that due to the finite hadronic radii the overlap of their wave functions should 
lead to effects of the Pauli principle on the quark level (cf. excluded volume) and to the appearance of 
multi-quark admixtures to the hadronic states [46 : |47| . When interpreted in terms of hadron-hadron 
scattering phase shifts (see, e.g., Refs. fisHiol: Is^" 3ll| for early examples of quark exchange contribu- 
tions) the quark substructure effects could be accounted for by a Beth-Uhlenbeck type generalization 
of the thermodynamical potential [sil; [s^; [s^- It is expected, that besides changes in the hadron wave 
functions, the masses of baryons shall get reduced at high densities whereas the masses of mesons 
could increase toward the chiral restoration. Such a behavior would distort the lines of constant chiral 
condensate towards the freeze-out data. Here, wc have disregarded the in-medium modifications of 
hadronic masses, which are currently being under consideration. 
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